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Abstract 

We study a model for leptogenesis in which the total CP asymmetries in the 
decays and scatterings involving the SU{2) singlet seesaw neutrinos N a vanish 
(e7v a = 0). Leptogenesis is possible due to non- vanishing CP violating lepton 
flavor asymmetries, realizing a situation in which the baryon asymmetry is 
due exclusively to flavor effects. We study the production of a net lepton 
asymmetry by solving the Boltzmann equations specific to this model, and 
we show that successful leptogenesis can be obtained at a scale as low as the 
TeV. We also discuss constraints on the model parameter space arising from 
current experimental upper limits on lepton flavor violating decays. 



1 Introduction 



Leptogenesis [TJ [2] (for a comprehensive review see ref . [3] ) is a theoretical mechanism 
that can explain the observed matter-antimatter asymmetry of the Universe. An 
initial lepton asymmetry, generated in the out-of-equilibrium decays of heavy singlet 
Majorana neutrinos, is partially converted in a baryon asymmetry by anomalous 
sphaleron interactions [I] that are standard model (SM) processes. Heavy Majorana 
singlet neutrinos are also a fundamental ingredient of the seesaw model [5], that 
provides an elegant explanation for the suppression of the neutrino masses with 
respect to all other SM mass scales. Leptogenesis can be quantitatively successful 
with a neutrino mass scale of the order of the atmospheric neutrino mass squared 
difference. This remarkable 'coincidence' links nicely the explanation of neutrino 
masses and of the baryon asymmetry within a single framework, and renders the 
idea that baryogenesis occurred through leptogenesis a very attractive one. 

In the standard seesaw case, computations of the CP violating (CPV) asymme- 
tries e7v Q in the decays of the N a singlet neutrinos include loop diagrams in which 
Majorana states appear in the internal lines [5J, and thus are lepton number vio- 
lating quantities [TJ. This is the reason why leptogenesis can proceed even when it 
is assumed that only one lepton flavor is relevant, as is the case at large tempera- 
tures (T > 10 12 GeV). However, at temperatures below ~ 10 12 GeV, lepton flavor 
dynamics plays an important role in leptogenesis, and cannot be neglected [SJ, E] 
(See [101 E] for earlier studies of flavor effects in leptogenesis, and [SJ d2] for recent 
reviews). In particular, in ref. [9] it was pointed out that leptogenesis can occur 
even when ctv q = 0, provided that the individual flavored CPV asymmetries e^-*^ 
(with j = e, /i, r) are non vanishing. 

Of course, e7v Q = Ylij e N a ^L 3 = means that total lepton number is not violated 
in N a decays. The reason why leptogenesis can still occur even in this case can be 
understood by analogy with the generation of a baryon asymmetry AB = B — B 
from a lepton asymmetry AL = L — L that, as is well known, does not require 
any baryon number violating CP asymmetry. Baryon number, or more precisely 
AB + AL, is in fact violated in the plasma by fast sphaleron reactions, with the 
result that part of AL is converted in AB yielding a ratio AB / AL = —28/51. 

Similarly, at T ~ Mjv a various interactions that are lepton and lepton flavor 
number violating occur in the plasma, like for example <-> N a <-> Of course, 
these reactions must be at least slightly out-of-equilibrium, otherwise they would 
quickly drive the individual ALj — > 0. However, the important point here is that 
generically these (washout) reactions proceed with different rates for different lepton 
flavors, erasing more efficiently, say, AL T than AL e41 . As a result, at T <C M Nl 
(Ni being the lightest heavy Majorana neutrino), once all washout processes are 
switched off, quite generically Ylj = 7^ results. This scenario, in which 
leptogenesis can proceed solely because of flavor effects, is what we call Purely 
Flavored Leptogenesis (PFL). It is worth noticing that PFL realizes the Sakharov 
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conditions [13] in a slightly different way that standard leptogenesis, since violation 
of lepton number occurs only in the washouts, while CP is violated only in the flavor 
charges, and the two conditions are thus disentangled. 

Ref. [H] analyzed the issue of the interplay between the lepton number breaking 
scale and the breaking scale of a flavor symmetry (of the Froggatt-Nielsen type [13]). 
It was found that in the case when the flavor symmetry is still unbroken during the 
leptogenesis era, but the vectorlike messengers masses are larger than the Majorana 
neutrino mass, the total CP asymmetry vanishes and a PFL scenario arises. In this 
paper we show that the PFL model of ref. [H] can indeed succeed in producing the 
cosmological baryon asymmetry. Interestingly, in this model there is an upper limit 
on the leptogenesis temperature fixed by the requirement that leptogenesis must 
occur in the flavored regime (T < 10 12 GeV) but, differently from the standard case, 
there is no lower limit and, as we will show, leptogenesis can be successful at a scale 
as low as the TeV. 

The rest of the paper is organized as follows: in section [2] we recall the main 
features of the model, we give the expressions for the flavor CPV asymmetries and 
we discuss an important rescaling property of the CPV asymmetries that leaves 
unaffected the washout rates. In section [3] we write down the Boltzmann Equations 
(BE) for the model and we present the main results. In section [5J we study some 
relations between the flavor violating parameters of the model and the low energy 
limits on lepton flavor violating processes. 

In our model washout and asymmetries in decays and scatterings occur at the 
same order in the couplings, and thus the derivation of the BE differs from the 
standard case in a non-trivial way. In appendix [A] we present a detailed derivation 
of the BE, that relies on the formalism introduced in ref. [16] to deal in a proper way 
with CPV asymmetries in scatterings. In appendix [B] we collect some definitions 
and useful formulae. 

2 The Model 

The model we consider here [H] is a simple extension of the SM containing a set 
of SU(2)l x 17(1) y fermion singlets, namely three right-handed neutrinos (N a = 
N a R + N^ R ) and three heavy vectorlike fields (F a = F aL + F aR ). In addition, we 
assume that at some high energy scale, taken to be of the order of the leptogenesis 
scale MjVj, an exact U(l)x horizontal symmetry forbids direct couplings of the 
lepton £i and Higgs $ doublets to the heavy Majorana neutrinos N a . At lower 
energies, U(l)x gets spontaneously broken by the vacuum expectation value a of 
a SU{2) L singlet scalar field S. Accordingly, the Yukawa interactions of the high 
energy Lagrangian read 

-C Y = l -N a M Na N a + F a M Fa F a + h ia lP R F a $> + N a (\ aa + \®>y B ) F a S + h.c. (1) 
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Figure 1: Effective seesaw operator for the light neutrino masses 



We use Greek indices a, (3 ■ ■ ■ = 1, 2, 3 to label the heavy Majorana neutrinos, Latin 
indices a,b - ■ ■ = 1, 2, 3 for the vectorlike messengers, and k, . . . for the lepton 
flavors e, /i, r. Following reference [H] we chose the simple U(l)x charge assignments 
X(£ Li , F La ,F Ra ) = +1, X(S) = -1 and X(N a , $) = 0. This assignment is sufficient 
to enforce the absence of iW$ terms, but clearly it does not constitute an attempt 
to reproduce the fermion mass pattern, and accordingly we will also avoid assigning 
specific charges to the right-handed leptons and quark fields that have no relevance 
for our analysis. The important point is that it is likely that any flavor symmetry 
(of the Froggatt-Nielsen type) will forbid the the same tree-level couplings, and 
will reproduce an overall model structure similar to the one we are assuming here. 
Therefore we believe that our results, that are focused on a new realization of the 
leptogenesis mechanism, can hint to a general possibility that could well occur also 
in a complete model of flavor. 

As it was discussed in [13], depending on the hierarchy between the relevant 
scales of the model (M Nl , M Fa , a), quite different scenarios for leptogenesis can 
arise. PFL arises when the relevant scales satisfy the hierarchy a < < Mp a 
that is, when the flavor symmetry U(l)x is still unbroken during the leptogenesis 
era and at the same time the messengers F a are too heavy to be produced in N\ 
decays and scatterings, and can be integrated away. As is explicitely shown by 
the last term in eq. ([I]), in general the vectorlike fields can couple to the heavy 
singlet neutrinos via scalar and pseudoscalar couplings, In ref. [2] it was assumed 
for simplicity a strong hierarchy A ^> A^ which allowed us to neglect all the A^. 
However, in all the relevant quantities (scatterings, CP asymmetries, light neutrino 
masses) at leading order the scalar and pseudoscalar couplings always appear in the 
combination A + A^ 5 **, and thus such an assumption is not necessary. The replacement 
A — ► A + A*- 5 - 1 would suffice to include in the analysis the effects of both type of 
interactions. 

2.1 Extended seesaw and light neutrino masses 

After U(l)x and electroweak symmetry breaking the Lagrangian eq. ([1]) generates 
masses for the light neutrinos through the effective mass operator depicted in fig- 



3 



ure [TJ The resulting mass matrix reads 
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where we have introduced effective seesaw-like couplings defined as 
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Note that, differently from standard seesaw, the neutrino mass matrix is of fourth 
order in the fundamental couplings (h and A) and includes an additional suppression 
factor of (a/M F ) 2 . 



2.2 N\ decays and CPV asymmetries 

Differently from standard leptogenesis in the present case, since Mp > M^, two- 
body Ni decays are kinematically forbidden. However, via off-shell exchange of the 
heavy F a fields, N\ can decay to the three body final states S$>1 and S<&1. The 
corresponding Feynman diagram is depicted in figure WLfl)- At leading order in 
r a = Mjsi 1 /Mp a) the total decay width reads [H] 

IV = J2T( Nl -> S% + m = ^ (^) 2 (AAt) n . (4) 
i ' 

As usual, CPV asymmetries in N\ decays arise from the interference between 
tree-level and one- loop amplitudes. As was noted in [Ti] , in this model at one- loop 
there are no contributions from vertex corrections, and the only contribution to 
the CPV asymmetries comes from the self-energy diagram W(b) ■ Summing over the 
leptons and vectorlike fields running in the loop, at leading order in r a the CPV 
asymmetry for N\ decays into leptons of flavor j can be written as 
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(5) 



Note that since the loop correction does not violate lepton number, the total CPV 
asymmetry that is obtained by summing over the flavor of the final state leptons 
vanishes j?], that is €\ = Y2j e ij — 0. This is the condition that defines PFL; namely 
there is no CPV and lepton number violating asymmetry, and the CPV lepton flavor 
asymmetries are the only seed of the Cosmological lepton and baryon asymmetries. 

It is important to note that the effective couplings A defined in eq. (J3]) are in- 
variant under the reparameterization 

A -> A • (rU)-\ -> (Ur) ■ , (6) 



4 



5 $ S $ 

! ! ! f>m ! 

U_L_ I r\ L 

$ 

(a) (6) 

Figure 2: Feynman diagrams responsible for the CPV asymmetry. 

where U is an arbitrary 3x3 non-singular matrix. Clearly the light neutrino mass 
matrix is invariant under this transformation. Moreover, also the flavor depen- 
dent washout processes, that correspond to tree level amplitudes that are deter- 
mined, to a good approximation, by the effective A couplings, are left essentially 
unchanged^] On the contrary, the flavor CPV asymmetries eq. (jSJ), that are deter- 
mined by loop amplitudes containing an additional factor of hr 2 h' , get rescaled as 
hr 2 h) — > h(rUr)' (rUr)fv . Clearly, this rescaling affects in the same way all the 
lepton flavors (as it should be to guarantee that the PFL conditions e a = e a j = 
are not spoiled), and thus for simplicity we will consider only rescaling by a global 
scalar factor r.U = U.r = re/ (with I the 3x3 identity matrix) that, for our pur- 
poses, is completely equivalent to the more general transformation (jH]). Thus, while 
rescaling the Yukawa couplings through 

A^Are" 1 , /i+ (7) 

does not affect neither low energy neutrino physics nor the washout processes, the 
CPV asymmetries get rescaled as: 

eij -> K 2 e Vj . (8) 

By choosing re > 1, all the CPV asymmetries get enhanced as re 2 and, being the 
Cosmological asymmetries generated through leptogenesis linear in the CPV asym- 
metries, the final result gets enhanced by the same factor. Therefore, for any given 
set of couplings, one can always find an appropriate rescaling such that the correct 
amount of Cosmological lepton asymmetry is generated. In practice, the rescaling 
factors re cannot be arbitrarily large: first, they should respect the condition that 
all the fundamental Yukawa couplings remain in the perturbative regime; second, 
as will be discussed in section [51 the size of the h couplings (and thus also of the 
rescaling parameter re) is also constrained by experimental limits on lepton flavor 
violating decays. 

1 The approximation is exact in the limit of pointlike F-propagators (s — Mp + iMpTp) — > Mp. 
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3 Boltzmann Equations 



In this section we compute the lepton asymmetry by solving the appropriate BE. In 
general, to consistently derive the evolution equation of the lepton asymmetry all 
the possible processes at a given order in the couplings have to be included. In the 
present case 1 <->• 3 decays and inverse decays, and 2 <-» 2 s, t and u channel scatter- 
ings all occur at the same order in the couplings and must be included altogether in 
the BE. The Feynman diagrams for these processes are shown in Figure [3j In ad- 
dition, the CPV asymmetries of some higher order multiparticle reactions involving 
the exchange of one off-shell Ni, also contribute to the source term of the asymme- 
tries at the same order in the couplings than the CPV asymmetries of decays and 
2^2 scatterings. More precisely, for a proper derivation of the BE it is essential 
that the CPV asymmetries of the off-shell 3 <-> 3 and 2 « 4 scattering processes 
depicted in figures and [8] in appendix [A] are also taken into account. In order 
to do this, we follow Ref. [16J and we split the BE for the evolution of the density 
asymmetry of the flavor ii as: 

Y ALi = (ValJi-s + (Y ALi h~2 + (W 3 t 3 + , (9) 

where Y AL . = (n^ — ng.)/s with (rig.) the number density of ii (anti)leptons, 
and s the entropy density. The time derivative is defined as Y = s H z dY j dz where 
z = M^/T and H is the Hubble parameter. The first term on the r.h.s. of eq. (Q 
represents the contribution of three body decays and inverse decays, the second term 
that of 2 <-> 2 scatterings, the third term is defined in terms of the off-shell (pole 
subtracted) 3^3 multiparticle density rates 7^3 = 73^3 — 73^ S 3 Cl1 , and similarly 
for the fourth term. The BE for Y At is derived by taking into account in full the first 
two terms on the r.h.s., while for the remaining two terms only the corresponding 
CP asymmetry is important, since non- resonant contributions to the washouts from 
multiparticle processes are always negligible. 

As regards the equation for the evolution of the heavy neutrino density , only 
the diagrams in fig. [3], that are of leading order in the couplings, are important. We 
refer to appendix [A] for a detailed derivation of the equations. The final result reads 

^ = -(^-1)7^ (10) 



Y A l, = (j/jvi - 1) e;7tot - Ay, 7* + (y Nl - 1) 7^ 



11) 



where in the last term of the_ second equation we have used the compact notation 
for the reaction densities 7^^ = j(Ni&i — > S&). ji and 7 tot are defined as 

7, = 7#,* + 7^ + 7^ + 7^ (12) 
7tot = Yl 7i + 7i, ( 13 ) 

i=e,[M,T 
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Figure 3: Feynman diagrams for 1 <-> 3 and 2 ^ 2 s, t and u channel processes after 
integrating out the heavy vectorlike fields F a . 

where in the second equation ji represents the sum of the CP conjugates of the 
processes summed in 7^. 

Since in this model Ni decays are of the same order in the couplings than scat- 
terings (that is C(A 2 )), the appropriate condition that defines the strong washout 
regime in the case at hand reads: 

— ^- > 1 (strong washout), (14) 
zHs z~i 

and conversely 7tot/(^-^ s )|«~i < 1 defines the weak washout regime. Note that this 
is different from standard leptogenesis, where at z ~ 1 two body decays generally 
dominate over scatterings, and e.g. the condition for the strong washout regime can 
be approximated as ^^/{zH s)\ z ~\ ~ Y^/ H\ z ^\ > 1. 



4 Results 

In this section we discuss a typical example of successful leptogenesis at the scale 
of a few TeV. The example presented is a general one. No particular choice of the 
parameters has been performed, except for the fact that the low energy neutrino 
data are reproduced within errors, and that the choice yields an interesting washout 
dynamics well suited to illustrate how PFL works. The numerical value of the final 
lepton asymmetry (Yal ~ —5.4 x 10~ 10 ) is about a factor of 3 larger than what is 
indicated by measurements of the Cosmic baryon asymmetry. This is however irrel- 
evant since, as was discussed in section 2., it would be sufficient a minor rescaling 
of the couplings (or a slight change in the CPV phases) to obtain the precise ex- 
perimental result. In the numerical analysis we have neglected the dynamics of the 
heavier singlet neutrinos since the N a masses are sufficiently hierarchical to ensure 
that A^2,3 related washouts do not interfere with iVi dynamics. Moreover, in the 
(strong washout) fully flavored regime (that is effective as long as T < 10 9 GeV) the 
iy 2j 3 CPV asymmetries do not contribute to the final lepton number asymmetry [18]. 

In figure H] we show the behavior of the various reaction densities for decays 
and scatterings, normalized to sHz, as a function of z = M^ 1 /T. The results 
correspond to a mass of the lightest singlet neutrino fixed to = 2.5 TeV, the 



7 




Figure 4: Reaction densities normalized to zHs for Ni — >• S£$ decays (red solid 
lines), s-channel SNi <-> scatterings (green dashed lines), and t,w-channel scat- 
terings in the point-like approximation (blue dotted lines). Left panel: r flavor. 
Right panel: electron flavor. 

heavier neutrino masses are M^ 2 = 10 TeV and Mn 3 = 15 TeV, and the relevant 
mass ratios r a = M^ x jMp a for the messenger fields are ri^s = 0.1,0.01,0.001 (the 
effects of the lightest F resonances can be seen in the s-channel rates in both panels 
in fig. H]). The fundamental Yukawa couplings h and A are chosen to satisfy the 
requirement that the seesaw formula eq. ([2]) reproduces within 2 a the low energy 
data on the neutrino mass squared differences and mixing angles [19J. Typically, 
when this requirement is fulfilled, one also ends up with a dynamics for all the 
lepton flavors in the strong washout regime. This is shown in the left panel in fig. [5] 
where we present the total rates for the three flavors. 

The left panel in fig. H] refers to the decay and scattering rates involving the 
T-flavor that, in our example, is the flavor more strongly coupled to Ni, and that 
thus suffers the strongest washout. It is worth noticing that, due to the fact that 
in this model scatterings are not suppressed by additional coupling constants with 
respect to the decays, the decay rate starts dominating the washouts only at z ^ 1. 
The right panel in fig. 0] depicts the reaction rates for the electron flavor, that is 
the more weakly coupled, and for which the strong washout condition eq. (j!4p is 
essentially ensured by sizeable s-channel scatterings. Scatterings and decay rates 
for the /i-flavor are not shown, but they are in between the ones of the previous two 
flavors. 

The total reaction densities that determine the washout rates for the different 
flavors are shown in the first panel in figure [51 The evolution of these rates with z 
should be confronted with the evolution of (the absolute value of) the asymmetry 
densities for each flavor, depicted in the second panel on the right. Since, as already 
stressed several times, PFL is defined by the condition that the sum of the flavor 
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CPV asymmetry vanishes (Y_V ey = 0), it is the hierarchy between these washout 
rates that in the end is the responsible for generating a net lepton number asymme- 
try. In the case at hand, the absolute values of the flavor CPV asymmetries satisfy 
the condition |e M | < |e e | < |e T |, as can be inferred directly by the fact that at z < 0.1, 
when the effects of the washouts are still negligible, the asymmetry densities satisfy 
this hierarchy. Moreover, since e Mje < while e r > 0, initially the total lepton num- 
ber asymmetry, that is dominated by Val t > is positive. As washout effects become 
important, the r-related reactions (blue dotted line in the left panel) start erasing 
Yal t more efficiently than what happens for the other two flavors, and thus the ini- 
tial positive asymmetry is driven towards zero, and eventually changes sign around 
z = 0.2. This change of sign corresponds to the steep valley in the absolute value 
I^alI that is drawn in the figure with a black solid line. Note that when all flavors 
are in the strong washout regime, as in the present case, the condition for the occur- 
rence of this 'sign inversions' is simply given by maxj gC]At ^|ej|/|Aij | 2 j ^ e r /|Ai r | 2 . 

From this point onwards, the asymmetry remains negative, and since the electron 
flavor is the one that suffers the weakest washout, 1al e ends up dominating all the 
other density asymmetries. In fact, as can be seen from the right panel in fig. it 
is Yal b that determines to a large extent the final value of the lepton asymmetry 
Y AL ="-5.4 x 10" 10 . 

A few comments are in order regarding the role played by the F a fields. Even 
if Mn 1 <C Mp a , at large temperatures z ^> 1 the tail of the thermal distributions 
of the N\, S and $ particles allows the on-shell production of the lightest F states. 
A possible asymmetry generated in the decays of the F fields can be ignored for 
two reasons: first because due to the rather large h and A couplings F decays 
occur to a good approximation in thermal equilibrium, ensuring that no sizeable 
asymmetry can be generated, and second because the strong washout dynamics 
that characterizes N\ leptogenesis at lower temperatures is in any case insensitive 
to changes in the initial conditions. 

In conclusion, it is clear from the results of this section that the model encounters 
no difficulties to allow for the possibility of generating the Cosmic baryon asymmetry 
at a scale of a few TeVs. Moreover, our analysis provides a concrete example of 
PFL, and shows that the condition e\ ^ is by no means required for successful 
leptogenesis. 

5 Lepton Flavor Violating Decays 

We have seen that a particular feature of this model is that the rescaling of the 
couplings eq. d7j) can enhance the CPV asymmetry by a factor k 2 (see eq. (jSJ)) 
without affecting neither the low energy neutrino physics nor the washout rates. 
In practice, it is this decoupling of the CPV asymmetries from the washouts that 
renders possible lowering the leptogenesis scale down to the TeV. It is then natural 
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Figure 5: Left panel: the total washout rates for each lepton flavor normalized 
to zHs as a function of z. Right panel: the evolution of the absolute value of 
the flavored density asymmetries and of the lepton number asymmetry (black solid 
line). The flavor CPV asymmetries are e le = —4.7 x 10~ 4 , = —1.9 x 10 -4 
and e lr = 6.6 x 10~ 4 . The final values of the asymmetry densities (at z > 1) are 
Y ALe = -7.1 x 1(T 10 , y AiM = -0.3 x 10- 10 , Y ALt = 2.0 x 10" 10 . 



to ask how large the rescaling factor k can be, or in other words how large the h 
couplings can become, without incurring in the violation of some phenomenological 
bound. To this aim, in this section we will derive upper bounds on the Yukawa 
couplings hi a from the non-observation of Lepton Flavor Violating (LFV) decays. 
The set of Yukawa interactions involving the heavy vectorlike fields F a , the SU(2) 
Higgs scalar <3> and lepton doublets can induce lepton flavor violating radiative 
decays £i — > £j~/. Here we will concentrate on /i — > cy that is the most constrained 
process. Note that, with respect to the radiative decays, lepton flavor violating 
decays like £' — > 3£ are more suppressed since they are induced by box diagrams 
rather then by penguin-type diagrams. Therefore we will not consider them. 
The partial decay width for the lepton flavor violating decay £i — > Iff reads 



T{£ t -> £ i7 ) 



a 



rri; 



10247T 4 Myy 



J2hlh ja F(M 2 w /M 2 Fa 



a=l 



(15) 



Here M w is the W ± SU{2) gauge boson mass which enter the loop through its 
longitudinal component $ ± and is the mass of the decaying lepton (the mass of 
the final state lepton rrij has been neglected). F(x) is a loop function given by 



F(x) 



x 



12(1 



x) 



(2 + 3x — 6a; 2 + x 3 + 6x log x) 



(16) 



and for x — > 0, F(x) — > x/6. Since Mp a S> Mw, F(M W /Mp a ) strongly suppresses 
the LFV radiative decays and thus, in general, the Yukawa couplings h ia will not 
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Figure 6: Br(fi — > ej) as a function of the mass of the lightest vectorlike field F±. 
No hierarchy for the h couplings and a mild hierarchy for the F masses (Mp 3 = 
1.3 • Mp 2 = 2 • Mf\) has been assumed. The region above the horizontal line is 
excluded by current experimental limits. The hatched region corresponds to non- 
perturbative couplings h > \f4n. 

be strongly constrained by the current experimental upper limits. In fact, as is 
shown in fig. [61 for Mp > a few TeVs and Yukawa couplings < .5, the radiative 
LFV decay rates always remain far below the present limits. Only for h ^ 2 (that 
is close to the perturbative limit represented by the hatched region in fig. [6]) and 
Mp 1 ~ 25 TeV (that was our choice in the numerical example) the fi — > ey decay 
rate becomes comparable to the present sensitivity, and is well within the reach of 
future experiments [20]. Thus we can safely conclude that present limits on LFV 
decays do not place any serious constraint on the viability of TeV scale leptogenesis 
within the PFL model discussed in this paper. 

6 Conclusions 

Variations of the standard leptogenesis scenario can arise from the presence of addi- 
tional (flavor) symmetries broken around or below the scale at which lepton number 
is effectively broken. Quite generically, the resulting scenarios can yield new qual- 
itative and quantitative changes on the way leptogenesis is realized. Here we have 
considered a model in which an Abelian U(l)x flavor symmetry, still unbroken dur- 
ing the leptogenesis era, is added to the SM gauge symmetry group. We have also 
assumed that the messengers fields responsible for the effective mass operators of the 
light particles are heavier that the lightest Majorana neutrino N± and thus, during 
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leptogenesis, cannot be produced. The model has the remarkable feature that the 
total CPV asymmetry in Aq decays vanishes, while the lepton flavor asymmetries 
are generically nonvanishing, and thus it constitutes an explicit realization of the 
scenario that we have called purely flavored leptogenesis. 

By using the BE specific for this model, we have studied the evolution of the 
asymmetry densities for the different lepton flavors as the temperature changes, 
and we have found that successful leptogenesis can occur at a scale as low as a 
few TeVs. This possibility is due to the fact that the size of the flavored CPV 
asymmetries is decoupled from the strength of the washouts and from low energy 
neutrino physics. This allows to rescale the CPV asymmetries up to rather large 
values, leaving unaffected the washout rates as well as the light neutrino masses and 
mixings. Our model shows that if new unbroken symmetries are present at a scale 
below the leptogenesis scale, this could have a very interesting and even surprising 
impact on the way leptogenesis is realized. 
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A Boltzmann Equations 

Following mainly ref. [16], we start by introducing a set of compact notations. We 
normalize particle densities to the equilibrium densities, y a = Y a /Y^ q where Y a = 
n a / s with n a the particle number density and s the entropy density, and we define 
the time derivative as Y = sHz dY j dz. Reaction densities are denoted by 7^ where 
A and B are respectively the initial and final states of the specific decay or scattering 
process. As in eq. (Q, we divide he BE for the evolution of the density asymmetry 
of the flavor ii Yal, = Ye i — Yg. into different contributions 

Y ALl = (Y ALi )i^3 + (Y ALi ) 2 ^ 2 + (Y ALi )^ + (Y ALi y 2 t, , (17) 
and we derive the explicit form of the different contributions in the following sections. 

A.l 1^3 and 2^2 processes 

The contributions (YalJi^>3 and (Y AL J 2 ^2 in e q- (ED arise from the 1 <-> 3 and 
2^2 reactions depicted in figure [31 that are all of 0(X 2 h 2 ). Using CPT invariance 
(1b = 1a) an d after linearizing in the CPV asymmetries and in the asymmetry 
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densities Ay a = y a — y a they can be written as: 

(^AiJiwsA = (Vm + 1)A 7 S* - (Ay 5 + Ay* + Ay07£*, (18) 

(*al« = (J/M + l)A 7 ^if - (j/^Aj/s + Ay* + A yii )j%l § , (19) 

(*AZ* W% = + l)A 7 ^f - (Ay s + y M Ay* + Ay 4 ) 7 Sf , (20) 

(^azJm^s = + 1)A75/ 1 - (Ays + Ay* + yNAydlS** , (21) 

where (Ya^W*-**^, (>azJ jv$^% and O'alJm^s are tne s ; * and u channel 
contributions to the 2^2 scattering term (Yal^*-^- For completeness, in these 
equations as well as in the following we keep trace of the density asymmetries Ay* 
and Ays of the Higgs and of the 5*-scalar. This is needed if one wishes to take into 
account spectator processes [21]. The related effects (that can be as large as 40% 
[2"T] ) depend, however, on the specific interactions of S with other particles (quarks) 
and are thus model dependent, and have been neglected in the present analysis. 
After summing up eqs. (I18H21I) we obtain 

(ValJi-s + (Y ALi ) 2 „ 2 = (y Nl + 1)(A 7 ^ + A 7 £f + A 7 £f + A<y£*) 

- Ay^ + t$ + >y%* + y Nll N s f) 

- A^( 7 S, + + y Nl iSl § + 7&*) 

" AjteCTffi. + y^Tg? + 7Sf + tS*) • (22) 

As is usual when only lowest order processes are included in the BE, the factor 
(y^Vi+1) signals an incorrect thermodynamical behavior (generation of an asymmetry 
in thermal equilibrium). In order to get the correct result we need to include in the 
BE also the CPV asymmetries of higher order processes, like the 3 «-> 3 and 2^4 
scatterings in which one Ni is exchanged as a virtual state in the internal lines. This 
is carried out in the following two sections. 



A. 2 3^3 processes 

Multiparticle tree level processes in which one Ni is exchanged in one internal line 
can be divided into on-shell and off-shell parts. For the on-shell parts, when the N\ 
line in the amplitude is cut, we obtain either the 1 «-> 3 or the 2^2 diagrams in 
fig. [3] which were already accounted for in eq. 022 1) . We then have to consider only 
the off-shell contributions, denoted as 7 ', where the superscript is a reminder that 
the given reaction has the on-shell piece subtracted out. The contribution of these 
off-shell processes to the washouts is negligible; however, the contribution of their 
CPV asymmetries cannot be neglected. Figure [7| shows the set of Feynman diagrams 
for 3 <-» 3 scattering processes yielding |ALj| = 2. They are of two types: ii 
and liii <-» 0. We do not show the analogous |ALj| = 1 diagrams, since they can 
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<I> 



s - 



■ S 



<i> 



■it 
-s 




Figure 7: Feynman diagrams for 3 3 scattering processes yielding |ALj| = 2. 
MiS «-> MiS (firstjour diagrams) ;_SjiS «-> (and <-> (fifth dia- 

gram); <-> f&iS'S' and £j£jS' <-> $$S* (last two diagrams). 



have either one £j or one £j (j ^ i) attached to one external leg, and for this reason 
their number is rather large. Since the CPV asymmetries of the full processes are of 
higher order in the couplings, at 0(X 2 h A ) the asymmetries of the off-shell parts are 
equal in magnitude and opposite in sign to the CPV asymmetries of their on-shell 
parts. In turn, the latter's are directly related to the CPV asymmetry in Ni decays 
and scatterings. Denoting the on shell parts of the rates as ^ on \ we then have: 

A7^ = -A 7 !T M (23) 
The on-shell pieces of the 3^3 reactions for |ALj| = 2 processes can be written as 

M«iS _ MiSpm , Mi p sNt , stu p *jvi , 5* pttN,. (2A) 

'<S>£iS ' Wi r MiS^ 'SNiQii ^ l$NiS£i ^ "liNi 5* ' 

(on) SiiS __ Sli pSNi (0 r\ 

'S£;S ' SNi tiS ' v ' 

The Pg represents the (decay or scattering) probability for the transition A — > £?, 
and were introduced in [16] to generalize the zero temperature branching ratios 
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to the finite temperature case, when the Nis have a finite probability to scatter 
inelastically with particles in the plasma before decaying. We have for example: 



Ni _ STVi 

PSl § = ^ P^ = ^ etc., (29) 

7tot 1 Ttot 

where 

7« = £ (iS* + iSU, + ifif + iff + + iff + iff + iff ) ■ (30) 

i 

Using the set of equations (12411280 together with the corresponding CP conjugate 
equations, including the contributions from |ALj| = 1 processes, and defining V£ = 
P 6 a + Pff, the different contributions to (^4)3^.3 can be written as: 



- A^v NlS V "P^ 5 - A^* 5 V P^ 

^7*4 ' $^ 1N1 ti 2_^t 1 *s 

(^a^Ws^^s = -Aj$j s p^ 1 

j 
j 

(^A^)i^s^i^5 = _ A7^. y^P^g 3 , 

J 

J 
j 



= -A7?E* - A ^ E^T 



w~«,* = - A 7 r E - A ^f E Kf , 

j j 

(iAXj^- .$«s*s = -A7 $€ . 2^ P *5 - MNiU2-j V *li ' 

(^w^ = - a 7 s* E ^ - A 7 j& E (si) 



15 



where in the l.h.s., of these equations, whenever the <-> £j transition is involved, it 
is left understood that % ^ j. To write the above relations, we have first completed 
the sums over flavor. For example for the second equation in f[3TT) : 



Ni 



{Y Mi U iS ^M jS = - A 7 ^ ! V sk + v Sk E A ^ 

= ~ A 7s<k E V sk + V sk E A ^ • (32) 



We have then used the specific relations valid for PFL (equivalent to ^ ■ e,- = 0) 

E ^sk = E A ^ = E A ^ = E A ^ = 0. (33) 
3 j j j 

Thus, the last term in the second line of eq. fl32l) vanishes, and the expression for 
(YM i )&t i s++&£ j s given in ([31]) is obtained. 

A. 3 2^4 processes 

The inclusion of 2 <-> 4 processes proceeds along similar lines than for 3^3 
processes. As is shown by the diagrams in figure [HI there are three types of |ALj| = 
2 contributions to the evolution equation for Yal^ The on-shell pieces of these 
contributions can be written as: 

(on)£ z S p i lS Ni (on) ii^ M Ni (on) $5 p J> S JVi (oa) 

By using again eq. ([231) to relate the on-shell parts to the relevant off-shell pieces, 
including the CP conjugate relations and including also the 2 <-> 4 processes with 
|ALj| = 1, we obtain the following contributions: 



w~*$ = - A ^f E ^ - A A E ^ 
(>alJ$4s#^5^ = -AifSiti E^S* ~ A ^S* E^k 

{YAL^-MiSS^li = - A lslti E^ 1 / ~~ A ^ 5 E 

C^)«^« OT = -At?^ £ ^ - A 7& E • (35) 



16 



Figure 8: Feynman diagrams for the 2^4 with |ALj| = 2 reactions S£i <-> <&S£i<&, 



i$ <-> SS^fc and $S «-> 



From the sets of equations (|3"I]) and ({35]), and using EjC^s + P£* + 

^'s* + ^s^ 3 ) = 1; we finally obtain the BE for the evolution of IatV 

Y ALi = (y Nl - 1)(A 7 ^ $ + A 7 £f + A 7 £f + At&*) (36) 

- Av.(t£* + 7 ff + vmiff + 7^) 

- Ay s (7£ # + y* l7 £ s + Tg,* + t£*) ■ 

where the term (y^ 1 — 1) in the r.h.s. shows that the correct thermodynamical 
behavior is recovered. 

Finally, using the approximate equalities between the scatterings and the decay 
asymmetries [16] 

A^i A/i* A<JViIi 

(37) 



A iVi A^ 1 * A'Y^ 15 ' An,* 1 ' 

TV! /V^ JViS iVil, 

7^*s 7^5 7^* 7$5 
the source term in equation (|36|) can be rewritten as 



A^ 1 

^al,) =(^-l)^P 7 ,, (3* 



if>S 



where 7i = 7 ^ 5 + 7 ^* + 7 ^ 5 + 7 $j^\ By using the relation 



N 1 



E^hs + liy Ei(7i + 7i)' 
and recalling that the flavored CPV asymmetries are given by 



(39) 



A 7 ? 



<i',S 



(40) 
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the evolution equation in (1381) can be recast as 



Y A Li =£i (ym - 1) 7tot - A?/j 7^ + (y Nl - 1) 7^ 



TV,* 



7* + (ym - 1)75/* 



Aye 



MS 



(41) 



where 7tot is given in eq. (1301) . Neglecting the terms in the second line proportional to 
Ay$ and A?/ s amounts to neglect 0(1) effects (21], and still yields a quite reasonable 
approximation. Finally, the evolution equation for the heavy Majorana neutrino can 
be simply written as: 

Y Nl = -(y Nl -l) ltot . (42) 



B Decay and Scatterings 

The thermally averaged reaction density for N\ <-» £i$>S is given by [21 ITT] 



Nl K,(z 



s ■ 



(43) 



where is the equilibrium number density for N±, Ki t2 {z) are Bessel functions 
(see appendix B in ref. [3]) and T^ s is the decay width given in PJ). The thermally 
averaged 2 — > 2 reaction densities are given by 



72^2 



M 4 Nl 



dxy/xKi {z\/~x) a 



512tt 5 2 J 1 

Here a(x) is the dimensionless reduced cross section 



(44) 



a(x) = ^r^* a ^ a ALA la F a (x) + 2 £ r a r fe Re[^^A; a A 16 G a ' 6 (x)]. (45) 



a<b 



The explicit form of the kinematical functions F a (x) and G a,b (x) depends on the 
specific s, t or u channel processes. In practice, for the t and u channels it is always 
a good approximation to use point-like interactions obtained by integrating out the 
heavy vectorlike fields. In this approximation, 7^* = 7$s^ = 72Z2 where 



t,u 

72U2 



M 4 Nl 
1024vr 5 2 



r a rb ~&Le(h* a hib\l a \ib) 



a<b 



/(*), (46) 



with 



/(*) 



dx 



x 



(47) 
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For s-channel scatterings (NiS <-> £i<&) the pointlike approximation is not sufficiently 
accurate, expecially at high temperatures T > M Nl , and the complete expression 
for the kinematical functions has to be used. We define 



x 



, x[(1-x)H$+(1 + x)H:$, 



with 



H$(x) 

H sS( X ) 



and 



2x(l - r 2 a x - 2ir a 7] a )(l - r\x + 2ir b r] b ), 

2r a r b + r a + r b + ir a r b (r] a - r] b ) + r a r b (x - 1), 

r a r b + r a + r b + ir a r b (r] a - r] b ) + (1 + ir a T] a )(l - ir b 7] b ) 

F?(x) = GT(x) . 



(48) 



(49) 
(50) 
(51) 

(52) 



The dimensionless parameter t] a introduced in the equations above corresponds 
to the total decay widths of the messenger fields F a normalized to the N\ neutrino 
mass: 



Va 



M Nl 



iirr n 



:i-r 2 a )(l + r a f\* la \ la + lj2 h *a h i 



(53) 
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